Rheology of Attractive Emulsions 
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We show how attractive interactions dramatically influence emulsion rheology. Unlike the re- 
pulsive case, attractive emulsions below random close packing, <j)RCP, can form soft gel- like elastic 
solids. However, above 4>RCP, attractive and repulsive emulsions have similar elasticities. Such com- 
pressed attractive emulsions undergo an additional shear-driven relaxation process during yielding. 
Our results suggest that attractive emulsions begin to yield at weak points through the breakage of 
bonds, and, above <j)RCP, also undergo droplet configurational rearrangements. 

PACS numbers: 83.60.-a, 83.80.Iz, 82.70.Kj, 61.43.-j 
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I. INTRODUCTION 



Emulsions are suspensions of droplets of one immis- 
cible fluid in another. They are widely used in tech- 
nological applications requiring the transport and flow 
of the dispersed fluid; these include oil recovery, food 
products, pharmacology, coatings, and cosmetics. The 
droplets are typically stabilized by a surfactant adsorbed 
on their interfaces; this ensures that there is a short- 
range repulsive interaction between them, which prevents 
their coalescence. Such a repulsive emulsion becomes a 
disordered elastic solid as it is compressed: this behav- 
ior can be characterized by the complex shear modulus, 
G*{uj) — G'{uj) + iG"(u)), where G' is the storage mod- 
ulus, G is the loss modulus, and U) is the angular fre- 
quency. For droplet volume fractions, 4>, approaching 
random close packing of spheres, (j)RCP ~ 0.64, from 
below, a repulsive emulsion exhibits a weak elasticity 
that arises from thermal fluctuations pj. By contrast, 
if <f> is increased above 4>rcpi the droplets are forced to 
deform; as a result, the elasticity is determined by the 
Laplace pressure scale of the droplets, a /a, where a is 
the interfacial tension between the dispersed and contin- 
uous phases and a is the average droplet size Even 
though it is an elastic solid, such an emulsion can nev- 
ertheless be made to flow quite easily: the imposition of 
sufficiently large shear causes it to yield. This behavior 
can be elucidated using oscillatory measurements of G' 
and G" performed at constant w and varying maximum 
strain amplitude 7. The viscoelastic moduli of a repul- 
sive emulsion above 4>rcp become strain-dependent for 
sufficiently large 7 as the emulsion yields. For increas- 
ing 7, the nonlinear G" exhibits a single, well-defined 
peak at a strain 7* before falling as ^~ u r ; Q' concomi- 
tantly decreases as j-"r with v' r ' « i/ r /2 [J fl H. The 
pronounced peak in G" is a characteristic feature of soft 



'Electronic address: weitz@seas.harvard.edu 



glassy materials; it is a direct consequence of a structural 
relaxation process and thus provides an effective way to 
characterize yielding (5!]. This approach is particularly 
useful for repulsive emulsions above 4>rcp, for which the 
peak in G" reflects the irreversible rearrangements of the 
densely-packed droplets [3, 0] . 

Another widely-encountered class of emulsions is char- 
acterized by droplets with additional attractive interac- 
tions between them. In stark contrast to the repulsive 
case, such an attractive emulsion can form an elastic 
solid even for <f> well below 4>rcp [2|; the bonds between 
droplets result in a connected network of aggregates that 
can support a shear stress @, H, OH [H|. As a result, 
an attractive emulsion must exhibit different flow and 
yielding behavior (l2| . However, despite its broad in- 
dustrial applications, exactly how this behavior occurs 
is unknown. Thus, measurements of the characteristic 7- 
dependent yielding of attractive emulsions are essential to 
elucidate how emulsion rheology depends on interdroplct 
interactions. 

In this Article, we explore the rheology of attractive 
emulsions using oscillatory measurements over a range 
of (f>. For increasing 7, G"{^f) of attractive emulsions 
below 4>rcp exhibits a single peak at a strain 7J <C 7* 
that increases with <f>. By contrast, G"{^) of attractive 
emulsions above 4>rcp exhibits two peaks at 7* and 72 ~ 
7*, unlike the repulsive case; these reflect two distinct 
structural relaxation processes. The time scales of both 
of these processes vary with shear rate as 7 _l/ with v k, 
0.8 — 1. Our results provide insight into the elasticity 
and yielding of attractive emulsions and highlight the 
sensitivity of emulsion rheology to attractive interactions. 



II. EXPERIMENTAL DETAILS 

We use emulsions comprised of silicone oil droplets 
dispersed in formamide, a solvent with negligible evapo- 
ration; the droplets are sterically stabilized by Pluronic 
PI 05, a non-ionic amphiphilic copolymer 13]. The me- 
chanical measurements are performed at T = 23° C on 
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strain- or stress-controlled rheometers (TA ARES G2 
or Anton Paar Physica MCR 501, respectively) using 
a parallel-plate geometry. The sample environment is 
controlled using a solvent trap. The plates are rough- 
ened to eliminate wall-slip, and we verify that measure- 
ments are independent of the gap size; furthermore, the 
results are similar to those obtained using a cone-and- 
plate geometry, indicating that they are not significantly 
influenced by the non-uniform strain field characteristic 
of a parallel-plate geometry. We pre-shear the samples 
prior to each measurement by imposing a 50s -1 constant 
shear-rate flow for 30s, followed by oscillatory shear of 
7 decreasing from 300% to 0.01% over a period of 230s. 
We do not observe sample creaming over the measure- 
ment duration. 

At the temperature and concentrations used here, 
P105 forms freely-dispersed globular micelles in for- 
mamide with radius a m — 6.5nm, aggregation number 
v m = 31 and critical micelle concentration c* = 20.2mM 
[l4| . These micelles induce attractive depletion interac- 
tions between emulsion droplets [3j , whose magnitude at 
interdroplet contact U we calculate using the Vrij model 
[l5| : U/UbT — Anaa^N 'a{c — c*)/v m , where a is the av- 
erage droplet radius, c is the P105 concentration, ks is 
Boltzmann's constant, T is temperature, and Na is Avo- 
gadro's number. These interactions preserve a lubricat- 
ing layer of formamide between the droplets, and hence 
are pairwise centrosymmetric and do not resist bend- 
ing. Because ions do not appreciably self-dissociate in 
formamide, and because we use a non-ionic surfactant, 
we do not expect electrostatics to play a significant role 

na. 

To investigate the influence of attraction on the rhe- 
ology of emulsions, we prepare repulsive emulsions with 
U = (c = 7mM < c*) and attractive emulsions with 
[7 = 7- 9k B T (c = 26.8 - 28.7mM) over a range of 
using high-shear rate homogenization or sonication with 
a probe tip. We densify the emulsions using centrifuga- 
tion and verify that they are stable against coalescence 
or Ostwald ripening (25[. We measure a = 6 — 9mN/m 
using a du Noiiy ring. The oil volume fraction <f> is deter- 
mined by measuring the masses of all components making 
up a sample. The droplets have a — 100, 106, or 128nm 
and relatively low polydispersity ~ 30 — 35% as mea- 
sured using dynamic light scattering. We also account 
for the thickness t of the surfactant layer adsorbed at the 
droplet surfaces to obtain the effective volume fraction 
(fi e ff « 0(1 +t/(a — t)) 3 . We estimate t 3.1nm based 
on neutron scattering data [TH, [TiJ and further verify this 
by fitting viscosity measurements of dilute repulsive sam- 
ples to simulation data appropriate for emulsions 0, Gil ■ 
This estimate for t is likely to be unchanged even at the 
highest 4> studied: the pressure required to compress the 
surfactant layer ~ ksT/s 3 , where its 2 is the interfacial 
area per polymer molecule, is approximately two orders 
of magnitude larger than the maximum osmotic pressure 
exerted on the droplets ~ O.la/a [H l2oj|. 

To test the generality of our results, we also study sili- 
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FIG. 1: Repulsive emulsions become fluid-like below 
random close packing. Linear storage and loss moduli, 
G'(oj) and G"(tu), of repulsive emulsions with U = 0, a = 
128nm, and $ e ff ~ 0.70, 0.65, and 0.60 (progressively lighter 
colors). 
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FIG. 2: Attractive emulsions are elastic both above 
and below random close packing. Linear storage and 
loss moduli, (a) G'(uj) and (b) G"(w), of attractive emulsions 
with U w 9k B T, a = lOOnm, and (j> eff w 0.73, 0.70, 0.66, 
and 0.62 (progressively lighter colors). For clarity, G' and 
G" data are multiplied by a factor of 0.5, 0.2, or 0.07 for 
<f> e ff = 0.70,0.66, or 0.62, respectively. 
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FIG. 3: Elasticities of attractive and repulsive emul- 
sions differ below random close packing but are sim- 
ilar above it. Linear storage and loss moduli, G' p (solid 
points) and G" (open points), measured at u> = 1 rad/s, 
normalized by the Laplace pressure a /a, for both attractive 
and repulsive emulsions of varying (j> e ff- Attractive emul- 
sions have a — 128nm, U ~ 9fcsT (circles), a — 106nm, 
U ~ IksT (upward triangles), and a — lOOnm, U ~ 9fcsT 
(squares); repulsive emulsions have a = 128nm (diamonds) 
and a = 106nm (downward triangles). The interfacial ten- 
sion a — 6 — 9mN/m for the range of surfactant concentra- 
tions used. Previous measurements of G' p (solid line) and G" 
(dashed line) for monodisperse repulsive emulsions [l], hori- 
zontally shifted to account for polydispersity, agree with our 
data. We find 4>rcp ~ 0.68 — 0.72, characteristic of 4>rcp for 
spheres of comparable polydispersity [171 ]. 
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FIG. 4: Repulsive emulsions above random close pack- 
ing exhibit one peak in G" '(7) at 7* while those below 
random close packing do not exhibit a peak in G"(7). 
Viscoelastic moduli G'(y) and G" (■y) measured at u — 1 rad/s 
of repulsive emulsions with U = 0, a = 128nm, and <f> e ff ~ 
0.70, 0.65, and 0.60 (progressively lighter colors). 



III. RESULTS AND DISCUSSION 



cone oil-in-water emulsions stabilized by an ionic surfac- 
tant, SDS. All of the data presented in this paper are for 
the oil-in-formamide emulsions stabilized by P105, with 
the exception of Figs. 7 and 8, which present data for 
the oil-in-water emulsions stabilized by SDS. These are 
well-defined and stable model systems characterized by 
very weak (< IksT) or strong (> lfc^T) attractive in- 
terdroplet interactions using an SDS concentration only 
slightly greater than or significantly greater than the crit- 
ical micelle concentration (c* ~ 8mM), respectively QHil- 
We prepare the emulsions using depletion fractionation; 
the droplets have a = 250nm and have low polydisper- 
sity ~ 10%. To investigate the influence of attraction on 
the rheology of emulsions, we study repulsive emulsions 
with U < lfcfiT (c = lOmM) and attractive emulsions 
with U w 2lk B T (c = 200mM). The effective volume 
fraction <j> e f f is defined to incorporate the effects of the 
thin liquid film between adjacent droplets due to their 
electrostatic repulsion. We study the emulsions using 
oscillatory rheology on a strain-controlled rheometer us- 
ing a roughened cone-and-plate geometry with a vapor 
trap. Further details and the rheology data for the SDS- 
stabilized repulsive emulsions are presented in [H, Hj] . 



The elastic behavior of emulsions is characterized by a 
w-independent regime of G'(u>); two examples, a repulsive 
emulsion with <f> e ff — 0.70 and an attractive emulsion 
with 4> e ff — 0.73, are shown in the top curves of Figs. 
1 and 2, respectively. Repulsive emulsions are solid-like 
when packed above 4>rcp but become fluid- like as 4> e ff 
is decreased below 4>rcp [Fig. 1] By contrast, at- 
tractive emulsions are solid-like over a wider range of 
4> e ff [Fig. 2]. To summarize this behavior, we plot the 
plateau modulus G', measured at lj = 1 rad/s, as a func- 
tion of 4> e ff- Above 4>rcp, G' p has the same magnitude 
for both attractive and repulsive emulsions as shown in 
Fig. 3. This indicates that the elasticity in both cases is 
dominated by the repulsive forces deforming the droplets. 
However, G" is an order of magnitude larger for attrac- 
tive emulsions, as compared to the repulsive case [Fig. 3]. 
For repulsive emulsions, G' p drops precipitously as 4> e ff is 
decreased below 4>rcp, indicating that the droplets are 
no longer compressed by contrast, the elasticity of 
attractive emulsions persists far below 4>rcp, as shown 
by the upper symbols in Fig. 3. 

To elucidate the microscopic mechanisms for emulsion 
flow, we investigate the yielding of emulsions with differ- 
ent <p e f f by measuring the 7-dependence of G' and G" at 
u = 1 rad/s. Repulsive emulsions above 4>RCP begin to 
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FIG. 5: Attractive emulsions above random close pack- 
ing exhibit two peaks in G"(7) at 7* and 72 while those 
below random close packing exhibit one peak at 7*. 

Viscoelastic moduli G'ip/) and G" (pf) measured at oj — 1 rad/s 
of attractive emulsions with U ~ 9fcsT and a — lOOnm (a- 
b) above 4>rcp (,4>&SS ~ 0.73 and 0.70, progressively lighter 
colors) and (c-d) below <j)RCP {4>elf ~ 0.66 and 0.62, progres- 
sively lighter colors). For clarity, G' and G" data are multi- 
plied by a factor of 0.7, 0.3, 0.8, or 0.2 for each 4>eff going 
from top to bottom. Straight lines indicate G' , G" ~ 7 _ " '~" 
for large 7. 




FIG. 6: Position of peaks in G"(7) for varying <f> e ff- 
Lower symbols represent 7*, upper symbols represent either 
72 or 7*- Attractive emulsions have a = 128nm, U ~ 9fcsT 
(circles), a = 106nm, U w 7ksT (upward triangles), and 
a — lOOnm, U « 9k bT (squares); repulsive emulsions have 
a = 128nm (diamonds) and a = 106nm (downward triangles). 



yield at 7* ~ 10% and G" exhibits a single, well-defined 
peak at this strain, as indicated by the arrow in Fig. 4 
[ll |J, ISj • Attractive emulsions above 4>rc p begin to yield 
at much smaller strain. Unexpectedly, G" exhibits two 
well-defined peaks, a first at ~ 1% <gC 7* and a sec- 
ond at 72 « 7*, before falling as 7"" , as indicated by 
the arrows in Fig. 5(b). Correspondingly, G' decreases 
weakly for 7 > 0.5% before falling as 7^ for 7 > 10%, 
with v" « v' /2 [Fig. 5(a)]. By contrast, G" of attractive 
emulsions below 4>rcp only exhibits a single peak at 7^, 
as shown in Fig. 5(d), and G' decreases smoothly for 
7 > 0.01-1% [Fig. 5(c)]. 

To understand these results, we examine the </> e //- 
dependence of 7*. Below 4>rcp, 7* decreases with de- 
creasing 4>eff, as shown by the lower symbols in Fig. 6 
[Lsl ] ; interestingly, this is similar to behavior predicted for 
particulate colloidal gels connected by "weak links" 0]. 
The elasticity of these emulsions results from the attrac- 
tive interdroplet bonds; these induce the formation of a 
stress-bearing connected network comprised of compact 
droplet aggregates @, [n], HH ■ We expect the elasticity 
to be dominated by the weakest bonds in the network, 
and hence yielding begins when these are broken. In this 
picture, a macroscopic deformation AL deforms such a 
bond by AL/(L/(), where L is the system size and ( is 
the characteristic distance between the weakest bonds. 
The force on such a bond is thus k a AL/(L/() = fc Q C7, 
where k a is the characteristic bond stiffness. Assuming 
that the bond breaks at a fixed critical force F* = fc a ^7^, 
the critical strain amplitude 7^ oc C • O ur measure- 
ments of 7^ thus imply that £ increases with decreasing 
4> e ff, suggesting that yielding begins at fewer, sparser 
weak points as <f> e f / is decreased. This may reflect the 
changing connectivity of the emulsion: the characteristic 
size £ of the droplet aggregates comprising the emulsion 
also increases with decreasing <f) e ff [j| [Tl|, [lj|. If the 
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FIG. 7: Volume fraction-dependent behavior is simi- 
lar for oil-in-water emulsions stabilized by SDS. (a) 

Linear plateau storage modulus, G' p 2], normalized by the 
Laplace pressure a/a and (b) position of peaks in G"(f), for 
varying <j> e ff. Emulsions have a = 250nm; attractive emul- 
sions have U ~ 21fcsT (circles) while repulsive emulsions have 
U < IksT (diamonds). The interfacial tension a = 9.8mN/m 
for the range of surfactant concentrations used. Upper sym- 
bols in (b) represent 7J while lower symbols represent 7*. 



FIG. 8: Strain-dependent behavior is similar for oil-in- 
water emulsions stabilized by SDS. Viscoelastic moduli 
(a) G'(-y) and (b) G"(7) of attractive emulsions stabilized by 
SDS with U ~ 21fc fl T and a = 250nm for <f> eff » 0.68, 0.65, 
0.62, and 0.59 (progressively lighter colors). Emulsions above 
4>RCP (top two curves) show peaks in G" at 71 and 72 while 
those below 4>rcp (lower two curves) show a single peak at 
7*. For clarity, G' and G" data are multiplied by a factor of 
0.7, 0.6, or 0.5 for <j) e ff = 0.65, 0.62, or 0.59, respectively. 



bonds between droplet aggregates are the weakest in the 
emulsion [20|, [H| , we would expect ( ps £ to increase with 
decreasing 4> e ff, consistent with our results. 

For increasing 7, attractive emulsions both below and 
above 4>rcp begin to yield and exhibit a peak in G" 
at 7 ~ 7^; in contrast, those above 4>rcp exhibit an 
additional peak in G" at 7 ~ 72 [Fig. 5]. This obser- 
vation implies the existence of an additional process by 
which the structure of these emulsions must relax before 
they can flow. This is likely due to the compression of 
the droplets: above (j)RCP, attractive emulsions require 
larger strain to fully yield through the irreversible re- 
arrangements of densely-packed droplets, just as in the 
repulsive case 0, H| . This hypothesis is supported by our 
observation that 72 ~ 7^ ■ The height of the first yielding 
peak in G"{"f) decreases as (p e ff increases above 4>rcp, 
as shown in Fig. 5(b), reflecting the increasing relative 
importance of these repulsive interactions as the droplets 
are increasingly compressed [22| . 

To test the generality of these results, we perform sim- 
ilar measurements on an oil-in-water emulsion stabilized 
by a different surfactant. The droplets are electrostat- 
ically stabilized by SDS, an ionic surfactant that also 



forms freely-dispersed globular micelles, and have low 
polydispersity (~ 10%) |l|,[3]. Similar to the case of P105- 
stabilized emulsions, G' p measured for the SDS-stabilized 
emulsions has the same magnitude for both attractive 
and repulsive emulsions packed above 4>rcp ~ 0.64 [Fig. 
7(a)]; this indicates that the elasticity is dominated by 
the repulsive forces deforming the droplets. For the re- 
pulsive emulsions, G' p drops precipitously as 4> e ff is de- 
creased below (j)RCP] by contrast, the elasticity of the 
attractive emulsions persists far below 4>rcp [Fig. 7(a)]. 
Moreover, G"(-f) of SDS-stabilized attractive emulsions 
above <pRCP exhibits two peaks at strains 7J ~ 1% < 7 r * 
and 72 ~ 7*, as indicated in Fig. 8(b). Correspondingly, 
G' decreases weakly for 7 > 0.5% before falling more 
quickly for 7 > 10% [Fig. 8(a)]. By contrast, G"(-y) 
of the attractive emulsions below 4>rcp exhibits a sin- 
gle peak at a strain 7*, and G' decreases smoothly for 
7 > 1% [Fig. 8]. These data are similar to those ob- 
tained for P105-stabilizcd emulsions; this confirms that 
our results are more general. 

The frequency-dependent mechanical response of 
emulsions directly reflects the time scales of their struc- 
tural relaxation. For repulsive emulsions above 4>rcp, 
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FIG. 9: (Color Online) Yielding of attractive emulsions 
is a shear-driven process. Constant 7 frequency sweep 
measurements of viscoelastic moduli G'(uj) and G"(uf) for 
two samples with a — lOOnm and U « 9}cbT above 4>RCP 
(a, 4>eff ~ 0.73) and below 4>rcp (b, 4>eff ~ 0.66), shifted 
onto a single master curve by normalizing by the shear-rate- 
dependent shift factors a(-y) and &(7). Inset: correspond- 
ing amplitude and frequency shift factors £1(7) (triangles) and 
6(7) versus 7 (squares). 



G"(oj) exhibits a shallow minimum at oj ~ 0.1 — 10 rad/s 
resulting from the combination of viscous loss at high uj 
and the configurational rearrangements of the individual 
droplets at low to [Fig. 1] [l|,l5j]. In contrast, for attrac- 
tive emulsions, G"(w) exhibits a minimum at uj m i n > 10 
rad/s, as shown in Fig. 2. This further suggests that 
attractive interactions alter the structural relaxation of 
emulsions. 

The structural relaxation process through which at- 
tractive emulsions yield occurs at frequencies much lower 
than is accessible in a linear rheological measurement, 
typical of soft glassy materials [||. We circumvent this 



problem by performing oj-dependent measurements hold- 
ing the strain rate amplitude 7 = ^uj constant for differ- 
ent values of 7. A relaxation process of time scale r 
gives rise to a peak in G" {10) at a frequency proportional 
to t _1 We observe one peak in G"(ui) for attractive 
emulsions with cp e ff < 4>rcp] strikingly, we observe two 
peaks in G"(ui) for those with <f> e ff > 4>RCP, as indi- 
cated in Fig. 9, in stark contrast to the repulsive case 
5]. This provides further confirmation that attractive 
emulsions above 4>r,cp undergo an additional relaxation 
process during yielding. 

The data measured at different 7 can be scaled onto a 
master curve by rescaling the moduli and frequency, as 
shown in Fig. 9. The frequency scaling factor 6(7) ~ r _1 
We explore its dependence on 7 to directly probe 
how the time scales of the structural relaxation processes 
of attractive emulsions depend on shear rate. We find 
6(7) ~ r _1 ~ 7" with v w 0.8 — 1 for both relaxation 
processes in attractive emulsions having 4> e f f both below 
and above <Prcp] two representative examples are shown 
in Fig. 9 (inset). This shear-driven behavior is similar 
to the yielding of colloidal gels |23f and to the yielding of 
repulsive emulsions above <Prcp [5]. 



IV. CONCLUSION 

Our data suggest a simple physical picture of emulsion 
flow. The rheology of emulsions with slippery bonds of 
magnitude U < 1&bT is similar to that of repulsive emul- 
sions; in stark contrast, emulsions with U > 7fcsT show 
a dramatically enhanced elasticity below 4>rcp- More- 
over, their nonlinear rheology is markedly different from 
the repulsive case. These attractive emulsions begin to 
yield under sufficient shear through the breakage of in- 
terdroplet bonds at weak points in the emulsion [T2I ]. 
Attractive emulsions above <Prcp also undergo shear- 
induced configurational rearrangements of the densely- 
packed droplets, similar to the repulsive case, at larger 
strain. This is reminiscent of the two-step yielding of at- 
tractive colloidal glasses . How this behavior changes 
with intermediate values of U remains to be explored. 

Emulsions are often used to model many other diverse 
physical systems and are crucial in various technologi- 
cal processes. We find that their linear and nonlinear 
rheology depend sensitively on the interactions between 
the droplets. Our results may thus be a useful input to 
such models and could help guide processes designed to 
control emulsion elasticity and flow behavior. 
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